Modelling of 3D fields due to ferritic inserts and test blanket modules in toroidal geometry at ITER

Computations in toroidal geometry are systematically performed for the plasma response to 3D magnetic perturbations produced by ferritic inserts (FIs) and test blanket modules (TBMs) for four ITER plasma scenarios: the 15 MA baseline, the 12.5 MA hybrid, the 9 MA steady state, and the 7.5 MA half-field helium plasma. Due to the broad toroidal spectrum of the FI and TBM fields, the plasma response for all the n = 1-6 field components are computed and compared. The plasma response is found to be weak for the high-n (n > 4) components. The response is not globally sensitive to the toroidal plasma flow speed, as long as the latter is not reduced by an order of magnitude. This is essentially due to the strong screening effect occurring at a finite flow, as predicted for ITER plasmas. The ITER error field correction coils (EFCC) are used to compensate the n = 1 field errors produced by FIs and TBMs for the baseline scenario for the purpose of avoiding mode locking. It is found that the middle row of the EFCC, with a suitable toroidal phase for the coil current, can provide the best correction of these field errors, according to various optimisation criteria. On the other hand, even without correction, it is predicted that these n = 1 field errors will not cause substantial flow damping for the 15 MA baseline scenario

Role of the pedestal position on the pedestal performance in AUG, JET-ILW and TCV and implications for ITER

The role of the pedestal position on the pedestal performance has been investigated in AUG, JET-ILW and TCV. When the pedestal is peeling-ballooning (PB) limited, the three machines show a similar behaviour. The outward shift of the pedestal density relative to the pedestal temperature can lead to the outward shift of the pedestal pressure which, in turns, reduces the PB stability, degrades the pedestal confinement and reduces the pedestal width. Once the experimental density position is considered, the EPED model is able to correctly predict the pedestal height. An estimate of the impact of the density position on a ITER baseline scenario shows that the maximum reduction in the pedestal height is 10% while the reduction in the fusion power is between 10% and 40% depending on the assumptions for the core transport model usedIn other plasmas, where the pedestal density is shifted even more outwards relative to the pedestal temperature, the pedestal does not seem PB limited and a different behaviour is observed. The outward shift of the density is still empirically correlated with the pedestal degradation but no change in the pressure position is observed and the PB model is not able to correctly predict the pedestal height. On the other hand, the outward shift of the density leads to a significant increase of eta(e) and eta(i) (where eta(e,i) is the ratio of density to temperature scale lengths, eta(e,i) = L-eta e,L-i/L-Te,L-i) which leads to the increase of the growth rate of microinstabilities (mainly ETG and ITG) by 50%. This suggests that, in these plasmas, the increase in the turbulent transport due to the outward shift of the density might play an important role in the decrease of the pedestal performance.Peer reviewe

A machine learning approach based on generative topographic mapping for disruption prevention and avoidance at JET

The need for predictive capabilities greater than 95% with very limited false alarms are demanding requirements for reliable disruption prediction systems in tokamaks such as JET or, in the near future, ITER. The prediction of an upcoming disruption must be provided sufficiently in advance in order to apply effective disruption avoidance or mitigation actions to prevent the machine from being damaged. In this paper, following the typical machine learning workflow, a generative topographic mapping (GTM) of the operational space of JET has been built using a set of disrupted and regularly terminated discharges. In order to build the predictive model, a suitable set of dimensionless, machine-independent, physics-based features have been synthesized, which make use of 1D plasma profile information, rather than simple zero-D time series. The use of such predicting features, together with the power of the GTM in fitting the model to the data, obtains, in an unsupervised way, a 2D map of the multi-dimensional parameter space of JET, where it is possible to identify a boundary separating the region free from disruption from the disruption region. In addition to helping in operational boundaries studies, the GTM map can also be used for disruption prediction exploiting the potential of the developed GTM toolbox to monitor the discharge dynamics. Following the trajectory of a discharge on the map throughout the different regions, an alarm is triggered depending on the disruption risk of these regions. The proposed approach to predict disruptions has been evaluated on a training and an independent test set and achieves very good performance with only one tardive detection and a limited number of false detections. The warning times are suitable for avoidance purposes and, more important, the detections are consistent with physical causes and mechanisms that destabilize the plasma leading to disruptions.Peer reviewe

Investigation of deuterium trapping and release in the JET divertor during the third ILW campaign using TDS

Selected set of samples from JET ITER-Like Wall (JET-ILW) divertor tiles exposed in 2015-2016 has been analysed using Thermal Desorption Spectrometry (TDS). The deuterium (D) amounts obtained with TDS were compared with Nuclear Reaction Analysis (NRA). The highest amount of D was found on the top part of inner divertor which has regions with the thickest deposited layers as for divertor tiles removed in 2014. This area resides deep in the scrape-off layer and plasma configurations for the second (ILW-2, 2013-2014) and the third (ILW-3, 2015-2016) JET-ILW campaigns were similar. Agreement between TDS and NRA is good on the apron of Tile 1 and on the upper vertical region whereas on the lower vertical region of Tile 1 the NRA results are clearly smaller than the TDS results. Inner divertor Tile 3 has somewhat less D than Tiles 0 and 1, and the D amount decreases towards the lower part of the tile. The D retention at the divertor inner and outer corner regions is not symmetric as there is more D retention poloidally at the inner than at the outer divertor corner. In most cases the TDS spectra for the ILW-3 samples are different from the corresponding ILW-2 spectra because HD and D-2 release occurs at higher temperatures than from the ILW-2 samples indicating that the low energy traps have been emptied during the plasma operations and that D is either in the energetically deep traps or located deeper in the sample.Peer reviewe

RF sheath modeling of experimentally observed plasma surface interactions with the JET ITER-Like Antenna

Waves in the Ion Cyclotron Range of Frequencies (ICRF) enhance local Plasma-Surface Interactions (PSI) near the wave launchers and magnetically-connected objects via Radio-Frequency (RF) sheath rectification. ITER will use 20MW of ICRF power over long pulses, questioning the long-term impact of RF-enhanced localized erosion on the lifetime of its Beryllium (Be) wall. Recent dedicated ICRF-heated L-mode discharges documented this process on JET for different types of ICRF antennas. Using visible spectroscopy in JET ICRF-heated L-mode discharges, poloidally-localized regions of enhanced (by similar to 2-4x) Be I and Be II light emission were observed on two outboard limiters magnetically connected to the bottom of the active ITER-Like Antenna (ILA). The observed RF-PSI induced by the ILA was qualitatively comparable to that induced by the JET standard, type-A2 antennas, for similar strap toroidal phasing and connection geometries. The Be II line emission was found more intense when powering the bottom half of the ILA rather than its top half. Conversely, more pronounced SOL density modifications were observed with only top array operation, on field lines connected to the top half of the ILA. So far the near-field modeling of the ILA with antenna code TOPICA (Torino Polytechnic Ion Cyclotron Antenna), using curved antenna model, was partially able to reproduce qualitatively the observed phenomena. A quantitative discrepancy persisted between the observed Be source amplification and the calculated, corresponding increases in E-// field at the magnetically connected locations to the ILA when changing from only top to only bottom half antenna operation. This paper revisits these current drive phased and half-ILA powered cases using for the new simulations flat model of the ILA and more realistic antenna feeding to calculate the E-// field maps with TOPICA code. Further, the Self-consistent Sheaths and Waves for Ion Cyclotron Heating Slow Wave (SSWICH-SW) code, which couples slow wave evanescence with DC Scrape-Off Layer (SOL) biasing, is used to estimate the poloidal distribution of rectified RF-sheath Direct Current (DC) potential V-DC in the private SOL between the ILA poloidal limiters. The approach so far was limited to correlating the observed, enhanced emission regions at the remote limiters to the antenna near-electric fields, as calculated by TOPICA. The present approach includes also a model for the rectification of these near-fields in the private SOL of the ILA. With the improved approach, when comparing only top and only bottom half antenna feeding, we obtained good qualitative correlation between all experimental measurements and the calculated local variations in the E-// field and V-DC potential.Peer reviewe

The effect of beryllium oxide on retention in JET ITER-like wall tiles

Preliminary results investigating the microstructure, bonding and effect of beryllium oxide formation on retention in the JET ITER-like wall beryllium tiles, are presented. The tiles have been investigated by several techniques: Scanning Electron Microscopy (SEM) equipped with Energy Dispersive X-ray (EDX), Transmission Electron microscopy (TEM) equipped with EDX and Electron Energy Loss Spectroscopy (EELS), Raman Spectroscopy and Thermal Desorption Spectroscopy (TDS). This paper focuses on results from melted materials of the dump plate tiles in JET. From our results and the literature, it is concluded, beryllium can form micron deep oxide islands contrary to the nanometric oxides predicted under vacuum conditions. The deepest oxides analyzed were up to 2-micron thicknesses. The beryllium Deuteroxide (BeOxDy) bond was found with Raman Spectroscopy. Application of EELS confirmed the oxide presence and stoichiometry. Literature suggests these oxides form at temperatures greater than 700 degrees C where self-diffusion of beryllium ions through the surface oxide layer can occur. Further oxidation is made possible between oxygen plasma impurities and the beryllium ions now present at the wall surface. Under Ultra High Vacuum (UHV) nanometric Beryllium oxide layers are formed and passivate at room temperature. After continual cyclic heating (to the point of melt formation) in the presence of oxygen impurities from the plasma, oxide growth to the levels seen experimentally (approximately two microns) is proposed. This retention mechanism is not considered to contribute dramatically to overall retention in JET, due to low levels of melt formation. However, this mechanism, thought the result of operation environment and melt formation, could be of wider concern to ITER, dependent on wall temperatures.Peer reviewe

Change in antihypertensive drug prescribing after guideline implementation: a controlled before and after study

<p>Abstract</p> <p>Background</p> <p>Antihypertensive drug choices and treatment levels are not in accordance with the existing guidelines. We aimed to assess the impact of a guideline implementation intervention on antihypertensive drug prescribing.</p> <p>Methods</p> <p>In this controlled before and after study, the effects of a multifaceted (education, audit and feedback, local care pathway) quality programme was evaluated. The intervention was carried out in a health centre between 2002 and 2003. From each health care unit (n = 31), a doctor-nurse pair was trained to act as peer facilitators in the intervention.</p> <p>All antihypertensive drugs prescribed by 25 facilitator general practitioners (intervention GPs) and 53 control GPs were retrieved from the nationwide Prescription Register for three-month periods in 2001 and 2003. The proportions of patients receiving specific antihypertensive drugs and multiple antihypertensive drugs were measured before and after the intervention for three subgroups of hypertension patients: hypertension only, with coronary heart disease, and with diabetes.</p> <p>Results</p> <p>In all subgroups, the use of multiple concurrent medications increased. For intervention patients with hypertension only, the odds ratio (OR) was 1.12 (95% CI 0.99, 1.25; p = 0.06) and for controls 1.13 (1.05, 1.21; p = 0.002). We observed no statistically significant differences in the change in the prescribing of specific antihypertensive agents between the intervention and control groups. The use of agents acting on the renin-angiotensin-aldosterone system increased in all subgroups (hypertension only intervention patients OR 1.19 (1.06, 1.34; p = 0.004) and controls OR 1.24 (1.15, 1.34; p < 0.0001).</p> <p>Conclusions</p> <p>A multifaceted guideline implementation intervention does not necessarily lead to significant changes in prescribing performance. Rigorous planning of the interventions and quality projects and their evaluation are essential.</p

Control of the hydrogen : deuterium isotope mixture using pellets in JET

Deuterium pellets are injected into an initially pure hydrogen H-mode plasma in order to control the hydrogen: deuterium (H:D) isotope mixture. The pellets are deposited in the outer 20% of the minor radius, similar to that expected in ITER, creating transiently hollow electron density profiles. A H: D isotope mixture of approximately 45%:55% is obtained in the core with a pellet fuelling throughput of Phi(pel) = 0.045P(aux)/T-e,T-ped similar to previous pellet fuelling experiments in pure deuterium. Evolution of the H: D mix in the core is reproduced using a simple model, although deuterium transport could be higher at the beginning of the pellet train compared with the flat-top phase.Peer reviewe

Current Research into Applications of Tomography for Fusion Diagnostics

Retrieving spatial distribution of plasma emissivity from line integrated measurements on tokamaks presents a challenging task due to ill-posedness of the tomography problem and limited number of the lines of sight. Modern methods of plasma tomography therefore implement a-priori information as well as constraints, in particular some form of penalisation of complexity. In this contribution, the current tomography methods under development (Tikhonov regularisation, Bayesian methods and neural networks) are briefly explained taking into account their potential for integration into the fusion reactor diagnostics. In particular, current development of the Minimum Fisher Regularisation method is exemplified with respect to real-time reconstruction capability, combination with spectral unfolding and other prospective tasks.Peer reviewe

Modification of the Alfven wave spectrum by pellet injection

Alfven eigenmodes driven by energetic particles are routinely observed in tokamak plasmas. These modes consist of poloidal harmonics of shear Alfven waves coupled by inhomogeneity in the magnetic field. Further coupling is introduced by 3D inhomogeneities in the ion density during the assimilation of injected pellets. This additional coupling modifies the Alfven continuum and discrete eigenmode spectrum. The frequencies of Alfven eigenmodes drop dramatically when a pellet is injected in JET. From these observations, information about the changes in the ion density caused by a pellet can be inferred. To use Alfven eigenmodes for MHD spectroscopy of pellet injected plasmas, the 3D MILD codes Stellgap and AE3D were generalised to incorporate 3D density profiles. A model for the expansion of the ionised pellet plasmoid along a magnetic field line was derived from the fluid equations. Thereby, the time evolution of the Alfven eigenfrequency is reproduced. By comparing the numerical frequency drop of a toroidal Alfven eigenmode (TAE) to experimental observations, the initial ion density of a cigar-shaped ablation region of length 4cm is estimated to be n(*) = 6.8 x 10(22) m(-3) at the TAE location (r/a approximate to 0.75). The frequency sweeping of an Alfven eigenmode ends when the ion density homogenises poloidally. Modelling suggests that the time for poloidal homogenisation of the ion density at the TAE position is tau(h) = 18 +/- 4 ms for inboard pellet injection, and tau(h) = 26 +/- 2 ms for outboard pellet injection. By reproducing the frequency evolution of the elliptical Alfven eigemnode (EAE), the initial ion density at the EAE location (r/a approximate to 0.9) can be estimated to be n(*) = 4.8 x 10(22) m(-3). Poloidal homogenisation of the ion density takes 2.7 times longer at the EAE location than at the TAE location for both inboard and outboard pellet injection.Peer reviewe